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Gluon propagator and confinement scenario in Coulomb gauge 
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We present numerical results in SU (2) lattice gauge theory for the instantaneous part of the gluon propagator 
in Coulomb gauge -D44,inst ~ K:oui(i?)<5(i)- Data are taken on lattice volumes 24* and 28* for 7 values of /3 in 
the interval 2.2 < /3 < 2.8. The data are confronted with the confinement scenario in Coulomb gauge. They are 
consistent with a linearly rising color-Coulomb potential Vcoui(i?). 



1. CONFINEMENT SCENARIO 
COULOMB GAUGE 



IN 



A particularly simple confinement scenario [y 
is available in the minimal Coulomb gauge. 
It attributes confinement of color to the long 
range of the color-Coulomb potential Vcoui(-R)- 
This quantity is the instantaneous part of 
the 4-4 component of the gluon propagator, 
D^^{x) = (gA^{x) gA^{0)), namely 1344 (x,t) = 
V^oui(|x|)(5(t) + P(x, t). The vacuum polarization 
term P(x, t) is less singular than d{t) at i = 0. 
Since A4 couples universally to color charge, the 
long range of T4oui {R) suffices to confine all color 
charge. It was conjectured that it is linearly ris- 
ing at large R, K;oui(^) ^ —o'couiR- If an exter- 
nal quark-antiquark pair is present, the physical 
potential Vw (R) between them may be extracted 
from a Wilson loop. The color-Coulomb potential 
contributes the term —CVcoui{R) directly to the 
Wilson loop, where C = {N"^ - 1)/{2N) in SU(N) 
gauge theory with external quarks in the funda- 
mental representation. The minus sign occurs be- 
cause the antiquark has opposite charge to the 
quark. The vacuum-polarization term is screen- 
ing, and one expects that Vw(^) is bounded 
above by this term asymptotically at large R, 
VwiR) < -CVcouiiR)- If ^w(^) is also finearly 
rising, Vw{R) ~ crR, where cr is the conventional 
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string tension, we get a < Cucoui- If dynamical 
quarks are present, the string "breaks" at some 
radius Rb, and the conventional asymptotic string 
tension vanishes, ct = 0. String-breaking is eas- 
ily explained in the Coulomb-gauge confinement 
scenario if K;oui(^^) is linearly rising even in the 
presence of dynamical quarks, as was also conjec- 
tured For, if so, it is energetically preferable 
to polarize a pair of sea quarks from the vacuum. 

Here we report the confrontation of this sce- 
nario with the numerical data for Vl:oui(^) for the 
case of pure gluodynamics. 

1.1. Relation of Koui to 

We may identify K;oui(-R) with the phenomeno- 
logical potential that is the starting point for 
QCD bound state calculations j^. The iden- 
tification of a phenomenological potential with 
the instantaneous part of the gluon propaga- 
tor, a fundamental quantity in the gauge the- 
ory, is possible because, remarkably, V^oui(-R) 
is a renormalization-group invariant, and thus 
scheme-independent, so it is independent of the 
cut-off A and of the renormalization mass /i. 
This follows from the non-renormalization of 
gAi, as expressed by the identity g(o)A_ 
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where and r refer to unrenormali- 



zed and renormalized quantities in the Coulomb 
gauge [0 . This identity has no direct analog in a 
Lorentz-covariant gauge. Because of the scheme- 
independence of Vcou\{R), its Fourier transform 
Vcoui(k) provides a scheme-independent defini- 
tion for the running coupling constant of QCD, 
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Table 1 



Values of fitting parameters from formula (|l|). 


a 




B 


A 




xVdof 


2.2 


7.38 ± 0.35 


8.15 ± 15.58 


8.83 ± 19.34 


9.91 ± 19.01 


0.81 


2.3 


5.81 ± 0.30 


24.57 ± 83.41 


32.57 ± 52.59 


17.26 ± 27.18 


0.72 


2.4 


6.60 ± 0.64 


6.72 ± 10.15 


7.70 ± 15.90 


3.99 ± 6.03 


0.20 


2.5 


7.45 ± 0.47 


6.00 ± 2.65 


6.16 ± 2.58 


3.06 ± 1.26 


0.05 


2.6 


8.23 ± 1.48 


8.61 ± 1.68 


4.53 ± 1.93 


3.84 ± 0.34 


0.07 


2.7 


11.30 ± 0.81 


7.73 ± 0.47 


6.87 ± 0.67 


3.64 ± 0.15 


0.07 


2.8 


8.18 ± 2.85 


10.38 ± 0.46 


0.58 ± 0.08 


0.06 ± 0.05 


0.19 



k2l/(k) = xo 5c'oui(|k|), and of a, ^ £!MA^. 
Here a;o — ii^^2N]^ ^"^^ ^com\ is a finite QCD 
mass scale ||]. 

2. NUMERICAL STUDY OF V^oui 

2.1. Method 

We have previously studied Q both space- and 
time-components of the gluon propagator, Dfj{k.) 

and Dii{k) = V^couKk) at equal time, at /3 = 2.2, 
on various lattice volumes, 14'* to 30"*, in the min- 
imal Coulomb gauge. This gauge is achieved nu- 
merically by (i) maximizing "Y^x i=i T^^^Ux,i with 
respect to all local gauge transformations g{x), 
where the sum is on all horizontal or spatial 
links, and then (ii) maximizing Ti:^Ux,4: with 
respect to all x-independent but a;4-dependent 
gauge transformations g{x4) where the sum is on 
all vertical or time-like links. This makes the 
3-vector potential Ai, for i = 1,2,3 transverse, 
^^A^ = 0, so A, = Af. 

In this gauge, the horizontal link variables Ux,i, 
for i = 1, 2, 3 are as close to the identity as possi- 
ble, but the vertical variables Ux,4 are much fur- 
ther from the identity. Not surprisingly, we found 
that, whereas £'*j(k) gave values that could be 
reasonably extrapolated to the continuum, this 
was not true for 1)44 (k). In order to remedy this, 
in the present study we extended our investiga- 
tion of £'44(k) to P = 2.2, 2.3, . . . , 2.8, on lattice 
volumes 24^ and 28**. We have also determined 
14oui(^) numerically by 2 quite different meth- 
ods. Method I relies on the standard formula 
i:»44(x,i) = Fcoui(|x|)5(0 -hF(x,i). Method II 



relies on the lattice analog of the continuum for- 
mula that is obtained |^ by integrating out A4 
namely, ^(|x-y|) = {{M-\~V^)M~^)\^^y). 
Here M{A) — —V • D(A) is the 3-dimensional 
Faddeev-Popov operator, and D(A) = V -|- Ax 
is the gauge-covariant derivative. Method II re- 
quires numerically inverting the lattice Faddeev- 
Popov matrix to obtain M~*(A)x,y, but it has 
the advantage that D44(k) is expressed entirely 
in terms of the spatial link variables Ux,i for 
i = 1,2,3, that are close to the identity. More- 
over, method II involves only the horizontal link 
variables that lie within a single time slice, so it 
is independent of the gauge fixing (ii) on verti- 
cal links. Thus it measures a truly instantaneous 
quantity. We found that method I did not exhibit 
scaling in the above range of /3, and we report here 
only the result of method II. 

2.2. Results 

For each /3 we have 100 configurations on vol- 
ume V = 24^ and 50 configurations on volume 
V — 28^. Runs have been done on the PC clus- 
ters at the Physics Department of New York Uni- 
versity and at the IFSC of Sao Paulo University. 

In figure 1 we plotted the results for f3 = 2.2, 
2.5, and 2.8 respectively. The horizontal axis 
measures k^ = 4a~^ sin^(n7r/L), for L = 24,28, 
rescaled to physical units by setting the physi- 
cal string tension equal to a — (0.44GeV)^ and 
using Table 3 of so the lattice spacing at 
P = 2.2, . . . , 2.8 is (in GeV"*) a = 1.066, 0.839, 
0.605, 0.433, 0.309, 0.231, 0.165, respectively. 
The vertical axis measures |k|'*V^oui(k) in phys- 
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ical units, so string tension may be read off from 
the vertical intercept (see below). Finite- volume 
artifacts are clearly visible at low momentum. To 
control these, in our fits we have dropped those 
low-momentum points for which appreciably dif- 
ferent values are obtained at volumes 24^ and 28* 
and, for the fit, we used only data points obtained 
&tV = 28*. 

A simple parametrization of Vcou\{R) would be 
—Vco\i\{R) — (^co\i\R — c/ R, which has the Fourier 
transform V;oui(k) = 87rcrcoui/|k|*-l-47rc/|k|2. We 
have used the fitting formulas 



IkrKoui(k) = A 

|k|*V;oui(k) = A 
|k|*V;oui(k) = A|k|^+Sk2 



Bk^ 



ln(l + k2/A2 



(1) 

(2) 
(3) 

Fit (|l|) has the asymptotic behavior at large k, 
consistent with the 1-loop /3-function. We re- 
port in Table |^ the values of the parameters from 
fit %. 

The most striking feature of the data is the fi- 
nite intercepts A. This is consistent with a finite 
string tension Ucoui = A/{?>tt). It scales rather 
nicely, varying from A — 7.38 ± 0.35 GeV^ at (3 = 
2.2 to A = 8.18 ± 2.85 GeV^ at /3 = 2.8, but with 
considerable variation in between. From the low- 
est and highest values, ^ = 5.8 and A — 11.3, we 
get respectively CTcouI = (0.48GeV)2 and ctcouI = 
(0.67GeV)2. The inequality a < (3/4)crcoui for 
SU(2), with cr = (0.44GeV)2, reads for these val- 
ues: 0.44 < 0.42(2) and 0.44 < 0.58(3). We ap- 
pear to be at or near saturation. 
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Figure 1. Fit of |k|*V;oui(k) using eq. (0). Data 
are for lattice volumes V — 24* ( x ) and V = 28* 
(*) and (3 = 2.2, 2.5 and 2.8 from top to bottom. 



